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Abstract: Cosmic ray antiprotons provide an important probe for the study of cosmic-ray propagation in tlie interstellar 
space and to investigate the existence of Galactic dark matter. The ARGO-YBI experiment, located at the Yangbajing 
Cosmic Ray Laboratory (Tibet, RR. China, 4300 m a.s.L, 606 g/cm^ ), is the only experiment exploiting the full coverage 
approach at very high altitude presently at work. The ARGO-YBJ experiment is particularly effective in measuring the 
cosmic ray antimatter content via the observation of the cosmic rays Moon shadowing effect. Based on all the data 
recorded during the period from July 2006 through November 2009 and a full Monte Carlo simulation, we searched for 
the existence of the shadow produced by antiprotons at the few-TeV energy region. No evidence of the existence of 
antiprotons was found in this energy region. Upper limits to the antip/p flux ratio are set to 5% at a median energy of 2 
TeV and 6% at 5 TeV with a confidence level of 90%. In the few-TeV energy range this result is the lowest available. 
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> 

o 
m 

o 



% 



1 Introduction 

Very High Energy Cosmic Ray (VHE CR) antiprotons are 
an essential diagnostic tool to approach the solution of 
several big topics of cosmology, astrophysics and parti- 
cle physics, besides for studying fundamental properties of 
the CR sources and propagation medium [ 1 ] . The enigma 
of the matter/anti-matter asymmetry in the local Universe, 
namely the existence of antimatter regions, the signatures 
of physics beyond the standard model of particles and 
fields, as well as the determination of the essential features 
of CR propagation in the insterstellar medium, are a few 
research topics which would greatly benefit from the de- 
tection of VHE antiprotons [2, 3, 4, 5, 6, 7, 8]. 

Cosmic rays are hampered by the Moon, therefore a deficit 
of CRs in its direction is expected (the so-called Moon 
shadow). Moreover, the Earth-Moon system acts as a mag- 
netic spectrometer. In fact, due to the Geomagnetic field 
(GMF) the Moon shadow shifts westward by an amount 
depending on the primary CR energy. The paths of primary 
antiprotons are therefore deflected in the opposite sense in 
their way to the Earth. This effect allows, in principle, the 
search for antiparticles in the opposite direction of the ob- 
served Moon shadow. 

If the energy is low enough and the angular resolution is 
good we can distinguish the two shadows, one shifted to- 
wards West due to the protons and the other shifted towards 
East due to the antiprotons. At high energies (> 10 TeV) 
the magnetic deflection is too small compared to the an- 



gular resolution and the two shadows cannot be disentan- 
gled. At low energies (<500 GeV) the shadows are much 
deflected but washed out by the poor angular resolution, 
thus the sensitivity is limited. Therefore, there is an opti- 
mal energy window for the measurement of the antiproton 
abundance. 

The ARGO-YBJ experiment is especially recommended 
for the measurement of the CR antimatter content via the 
observation of the Galactic CR shadowing effect due to: 
(1) good angular resolution and pointing accuracy and their 
long term stability; (2) low energy threshold; (3) real sen- 
sitivity to the GMF. Indeed, the low energy threshold of 
the detector allows the observation of the shadowing effect 
with a sensitivity of about 9 standard deviations (s.d.) per 
month at TeV energy. 

In this paper we report the measurement of the CR p/p flux 
ratio in the TeV energy region with all the data collected 
during the period from July 2006 to November 2009. 

2 The ARGO-YBJ experiment 

The detector is composed of a central carpet large ^74 x 
78 m^, made of a single layer of Resistive Plate Chambers 
(RPCs) with ^-^93% of active area, enclosed by a guard ring 
partiaUy (^20%) instrumented up to ^lOOx 1 10 vc?. The 
apparatus has modular structure, the basic data acquisition 
element being a cluster (5.7x7.6 m^), made of 12 RPCs 
(2.8 X 1.25 m^ each). Each chamber is read by 80 external 
strips of 6.75x61.8 cm^ (the spatial pixels), logically or- 
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ganized in 10 independent pads of 55.6x61.8 cvc? which 
represent the time pixels of the detector The read-out of 
18360 pads and 146880 strips are the experimental output 
of the detector [9]. The central carpet contains 130 clusters 
(hereafter, ARGO-130) and the full detector is composed 
of 153 clusters for a total active surface of '^6700 vc? . 

The whole system, in smooth data taking since July 2006 
firstly with ARGO-130, is in stable data taking with the 
full apparatus of 153 clusters since November 2007 with 
the trigger condition Nfrig = 20 and a duty cycle >85%. 
The trigger rate is ^3.5 kHz with a dead time of 4%. 

The analysis reported in this paper refers to events collected 
after the following selections: (1) more than 20 strips on 
the ARGO-130 carpet; (2) zenith angle of the shower ar- 
rival direction less than 50°; (3) reconstructed core posi- 
tion inside a 150x150 m^ area centered on the detector 
According to the simulation, the median energy of the se- 
lected protons is E50 «1.8 TeV (mode energy «0.7 TeV). 

3 Monte Carlo simulation 

A detailed Monte Carlo (MC) simulation of the CR propa- 
gation in the Earth-Moon system [10] has been developed 
in order to estimate the expected antiproton flux in the side 
opposite to the CR Moon shadow. The simulation is based 
on the real data acquisition time. The Moon position has 
been computed at fixed times, starting from July 2006 up 
to November 2009. Such instants are distant 30 seconds 
each other. For each time, after checking that the data ac- 
quisition was effectively running and the Moon was in the 
field of view, primaries were generated with arrival direc- 
tions sampled within the Moon disc basing on the effective 
exposure time. 

After accounting for the arrival direction correction due to 
the magnetic bending effect, the air shower development 
in the atmosphere has been generated by the CORSIKA 
V. 6.500 code with QGSJET/GHEISHA models [11]. CR 
spectra have been simulated in the energy range from 10 
GeV to 1 PeV following the results given in [12]. About 
10^ showers have been sampled in the zenith angle interval 
0-60 degrees. The experimental conditions (trigger logic, 
time resolution, electronic noises, etc.) have been repro- 
duced by a GEANT4-based code [13]. 

4 Data analysis 

For the analysis of the shadowing effect, the signal is col- 
lected within a 10° x 10° sky region centered on the Moon 
position. We used celestial coordinates (right ascension 
and declination, hereafter R.A. and DEC.) to produce the 
event and background sky maps, with 0.1°x0.1° bin size. 
Finally, after a smoothing procedure, the significance map, 
used to estimate the statistical significance of the observa- 
tion, is built. 





xlQ-* 










684 


— 










682 
680 


'- 


I. 


iti 


pBm'J r l™lEt 


i 


678 


'- 


m 


T 




mi 


676 


- J 


w ^ 




Wr 


674 


'-m 


1 






Tfc 


672 










fM 


670 








■ Events 


il. 


668 


- 






■ Background 




,,,,,,, Till 



-10 -8 -6 -4 -2 2 4 6 8 10 
w E 

Figure 1 : Deficit of CRs around the Moon position pro- 
jected along the R.A. direction. Showers with N>60 
recorded from July 2006 until November 2009 are shown. 

As can be noticed from Fig. 1, the Moon shadow turns out 
to be a lack in the smooth CR signal observed by ARGO- 
YBJ, even without subtracting the background contribu- 
tion. The background events are not uniformly distributed 
around the Moon, because of the non-uniform exposure of 
the map bins to CR radiation. 

The background has been estimated with two different 
approaches, the time swapping and the equi-zenith angle 
methods, as described in [14], in order to investigate possi- 
ble systematic uncertainties in the background calculation. 
They give results consistent with each other within 1 s.d.. 
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Figure 2: Moon shadow significance map. It collects all 
the events detected by ARGO-YBJ from July 2006 until 
November 2009. The event multiplicity is 20<N<40 and 
zenith angle < 50°. The color scale gives the statistical 
significance. 

A significance map of the Moon region is shown in Fig. 2. 
It contains all the events belonging to the lowest multiplic- 
ity bin investigated (20<N<40), collected by ARGO-YBJ 
during the period July 2006 - November 2009 (about 3200 
hours on-source in total). The significance of the maximum 
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Figure 3: Deficit counts measured around the Moon pro- 
jected along the East-West axis for different multipHc- 
ity bins (black circles) compared to MC expectations (red 
squares). 

deficit is about 22 s.d.. The observed westward displace- 
ment of the Moon shadow by about 1.5° allows to appre- 
ciate the sensitivity of the ARGO-YBJ experiment to the 
GMF. In such a map a potential antiproton signal is ex- 
pected eastward within 1.5° from the actual Moon position 
(i.e., within 3° from the observed Moon shadow). The me- 
dian energy of the selected events is E50 «750 GeV (mode 
energy w 550 GeV) for proton-induced showers. The cor- 
responding angular resolution is ^1.9°. However, the large 
displacement of the shadow is only one ingredient of this 
analysis, the other being the angular resolution which is not 
small enough in this multiplicity range. Indeed, as can be 
seen from Fig. 2, the matter shadow extends to the antimat- 
ter side with a significance of about 10 s.d., thus limiting 
the sensitivity to the antiproton abundance measurement. 

We stress that this is the first time that an HAS array is able 
to detect the Moon shadow so shifted, observing the signal 
due to sub-TeV primary CRs. 

5 Results and discussion 

The optimal energy windows for the measurement of the 
antiproton abundance in CRs are given by the following 
multipHcity ranges: 40<N< 100 and N> 100. In the former 
bin the statistical significance of Moon shadow observation 
is 34 s.d., the measured angular resolution is ^1°, the pro- 
ton median energy is 1.4 TeV and the number of deficit 
events about 183000. The Moon shadow is shifted by 0.82° 
westward from its actual position. In the latter multiplicity 
bin the significance is 55 s.d., the measured angular reso- 
lution ^0.6°, the proton median energy is 5 TeV and the 
number of deficit events about 46500. The Moon shadow 
is shifted by 0.30° westward. 

The chance of unfolding all CR spectral contributions relies 
on MC simulations, as well as the search for antiprotons 
demands to properly reproduce the Moon shadow signal. 
The deficit counts observed around the Moon projected on 



the East- West axis are shown in Fig. 3 for 4 multiplicity 
bands compared to MC expectations. We used the events 
contained in an angular band parallel to the East- West axis 
and centered on the observed Moon position. The widths of 
these bands are chosen on the basis of the MC simulation 
so that the shadow deficit is maximized. They turn out to 
be proportional to the Nst^ip-dependent angular resolution: 
±2.9° in 40< Nst„p < 60, ±2.6° in 60< Nstrip < 100, 
±2.1° in 100< N,trip < 200, ±1.6° in 200< N^trip < 
500. As an expected effect of the GMF, the profile of the 
shadow is broadened and the peak positions shifted west- 
ward as the multiplicity (i.e., the cosmic ray primary en- 
ergy) decreases. The data are in good agreement with the 
expectations from the MC simulation. 

In order to evaluate the p/p ratio, the projections of the 
Moon shadow along the R.A. direction have been consid- 
ered. The GMF shifts westward the dip of the signal from 
positively charged primaries. Searching antiprotons means 
looking for excesses in the eastern part of the R.A. pro- 
jection, i.e. trying to fit the Moon shape expected from 
combining CR and antimatter to the shape obtained from 
the experimental data. It is worth noticing that for the an- 
tiprotons we assume the same energy spectrum of the pro- 
tons. Of course, whichever matter-antiprotons combination 
is obtained, the total amount of triggered events must not 
be changed, so that the fitting procedure consists in trans- 
ferring MC events from the CRs to the antiproton shadow 
and comparing the result with the data. 

To make such a comparison, we firstly adopted the follow- 
ing method. We obtained two kinds of Moon shadow, cast 
by all CRs and protons, respectively. After projecting them 
along the R.A. direction, we used a superposition of several 
Gaussian functions to describe the deficit event distribution 
in each shadow. Four Gaussian functions were found to be 
adequate for fitting both distributions within 5° from the 
Moon disc center Let us name 9 the angular distance from 
the Moon disc center and fm{9) the Gaussian function su- 
perposition describing the CR shadow. Let J-p{9) be the 
proton shadow, obtained by imposing a given power law 
spectrum. The observed Moon shadow should be expressed 
by the following function: 

hiooNie) = {l-r) f,n{0) +rTp{e) 
= {l-r)f„,i9)+rTp{-9) 

(0 < r < 1) where the first term represents the deficit in 
CRs and the second term represents the deficit in antipro- 
tons. This function must be fitted to the data to obtain the 
best value of r. 

We also applied a second method to determine the an- 
tiproton content in the cosmic radiation. Without introduc- 
ing functions to parameterize the expectations, we directly 
compared the MC signal with the data. We performed a 
Maximum Likelihood fit using the p content as a free pa- 
rameter with the following procedure: 

1. the Moon shadow R.A. projection has been drawn 
both for data and MC. 
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2. the MC Moon shadow has been split into a "matter" 
part plus an "antiproton" part, again so that the total 
amount of triggered events remains unchanged: 

^Mc{mat) — > ^Mc{r;mat + p) = 
= (1 - r)^Mc{mat) + ^mc{p) 

3. for each antiproton to matter ratio, the expected 
Moon shadow R.A. projection ^Mci''"-, fncit + p) is 
compared with the experimental one via the calcula- 
tion of the likelihood function; 

B 

logC{r) = J2 NMEiir)] - E,{r) - ln{K^:) 

4=1 

where N^ is the number of experimental events in- 
cluded within the i-th bin, while Ei (r) is the num- 
ber of events expected within the same bin, which 
is calculated by adding the contribution expected 
from MC {^Mc{f\ mat+p)) to the measured back- 
ground. 

Both methods described above give results consistent 
within 10%. The r parameter which best fits the expec- 
tations to the data turns out to be always negative, i.e. it 
assumes non-physical values throughout the whole energy 
range investigated. With a direct comparison of the R.A. 
projections, the r-values which maximize the likelihood 
are: -0.076±0.040 and -0.144±0.085 for 40<N<100 and 
N>100, respectively. The corresponding upper limits with 
90% confidence level (c.l.), according to the unified Feld- 
man & Cousins approach [15], are 0.034 and 0.041, respec- 
tively. Since the anti-shadow was assumed to be the mirror 
image of the proton shadow, we assume for the antipro- 
tons the same median energy. The p/p ratio is ^{p)/^{p) 
= l/fp- $(p)/$(matter), therefore, being the assumed 
proton fraction /p=73% for 40<N<100 and /p=71% for 
N>100 [12], we obtain the following upper limits at 90% 
C.L: 0.05 for 40<N<100 and 0.06 for N> 100. Notice that 
the two values are similar, in spite of the different multi- 
plicity interval. It is a consequence of the combination of 
the two opposing effects of the angular resolution and of 
the geomagnetic deviation. 

In Fig. 4 the ARGO-YBJ results are shown with all the 
available measurements. The solid curves refer to a theoret- 
ical predictions for a pure secondary production of antipro- 
tons during the CR propagation in the Galaxy by Donato 
et al. [7]. The curves was obtained using the appropriate 
solar modulation parameter for the PAMELA data taking 
period [5]. The long-dashed lines refer to a model of pri- 
mary p production by antistars [8]. The rigidity-dependent 
confinement of CRs in the Galaxy is assumed oc R^*, and 
the two curves represent the cases of S = 0.6, 0.7. The dot- 
dashed line refers to the contribution of p from the annihi- 
lation of a heavy dark matter particle [6]. The short-dashed 
line shows the calculation by Blasi and Serpico [16] for 
secondary antiprotons including an additional p component 
produced and accelerated at CR sources. 
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Figure 4: The p/p flux ratio obtained with the ARGO-YBJ 
experiment compared with all the available measurements 
and some theoretical predictions (see text). 

6 Conclusions 

The ARGO-YBJ experiment is observing the Moon 
shadow with high statistical significance at an energy 
threshold of a few hundred GeV. Using all data collected 
until November 2009, we set two upper limits on the p/p 
flux ratio: 5% at an energy of 1.4 TeV and 6% at 5 TeV 
with a confidence level of 90%. In the few-TeV range the 
ARGO-YBJ results are the lowest available, useful to con- 
strain models for antiproton production in antimatter do- 
mains. 



References 

[1] F.W. Stecker, in "Matter-Antimatter Asymme- 
try" ed. L. Iconomidou-Fayard & J. Tran Thanh 
Van (Hanoi:The Gioi) pp. 5-14 (2003) (preprint 
arXiv:hep-ph/0207323v 1 ) 
[2] FW. Stecker and Wolfendale A.W., Nature 309 

(1984) 37. 
[3] M.Yu. Khlopov et al., Phys. Rev. D62 (2000) 083505. 
[4] A.W. Strong et al., Annu. Rev Nucl. Part. Sci. 57 

(2007) 285. 
[5] O. Adriani et al., Phys. Rev Lett. 105 (2010) 121101. 
[6] M. Cirelli et al., Nucl. Phys. B813 (2009) 1. 
[7] F. Donato et al., Phys. Rev Lett. 102 (2009) 071301. 
[8] S.A. Stephen and R.L. Golden, Space Sci. Rev. 46 

(1987)31. 
[9] G. Aielh et al., NIM A562 (2006) 92. 
[10] G. Di Sciascio and R. luppa, NIM A630, 301 (201 1). 
[11] D. Heck et al.. Report FZKA 6019 (1998). 
[12] J.R. Horandel, Aifrop. Phys 19 193, 2003. 
[13] YQ. Guo et s\.,Chinese Phys. C34 (2010) 555. 
[14] G. AielH et al., ApJ, lU L208, 2010. 
[15] G. Feldman and R. Cousins, Phys. Rev D57 (1998) 

3873. 
[16] R Blasi and RD. Serpico Phys. Rev Lett. 103 (2009) 
081103. 



